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1. Introduction 

Switching on the coffee machine, how long does it take until you smell the coffee for 
the first time? In statistical physics questions of this type are known as first-passage 
problems [1]. Equivalently one may ask; for the persistence probability that someone, 
after turning on the coffee machine, does not smell the coffee for a certain period of 
time. In more general terms, the persistence probability is defined as the probability 
that a certain observable in a random process does not cross its expection value from 
the initial state until time t. 

In the last decade persistence has been studied both theoretically and 
experimentally in a large variety of systems [2-75]. There are two main categories of 
persistence studies. For global observables such as the total magnetization or the 
average particle density, one studies the global persistence probability Pg{t) that an order 
parameter does not cross its expectation value up to time t [11-26]. In critical dynamical 
systems it is observed that the global persistence probability decays algebraically as 
Pg{t) ~ t'^^. For Markovian dynamics one can show that the global persistence 
exponent can be expressed by a scaling relation 6gZ — \ — d -\- 1 -\- r]/2, where d is 
the spatial dimension while 77 and A are the static and the autocorrelation exponents, 
respectively. 

The other category of persistence studies, on which we will focus in the present 
work, deals with local observables such as the occupation number or spin orientation 
at individual sites. Here, one is interested in the local persistence probability Pe{f,t) 
that the local state at a given position r in space has not changed up to time t. For 
example, in spin models Pe{f, t) may be defined as the probability that a spin at a given 
site does not flip until time t. If the system is translationally invariant, the persistence 
probability does not depend on r so that Ptit) is just the fraction of spins that did not 
flip until time t. 

In many systems evolving towards a scale-free state the local persistence probability 
is found to decay algebraically as 

Pt{t) ~ t-'^ , (1) 

where is the so-called local persistence exponent [2-8,21,23,27-35]. This exponent is 
generally different from 9g and seems to have certain universal features, i.e., its numerical 
value may coincide in various models with continuous phase transitions belonging to the 
same universality class. However, in contrast to 6g, which can be expressed in terms of 
the bulk exponents by a scaling relation, 9(, seems to be an independent exponent. So 
far there are only very few exact results [3,5-8,36,37]. One of the most pathbreaking 
results was obtained by Derrida et al. [6,7], who were able to compute the persistence 
exponent 9^ in the one-dimensional Ising and Potts model with heat bath dynamics. 
This exact solution maps the Glauber-Ising dynamics to a time-reversed coagulation 
process with a special boundary condition and relates 6*^ to a universal amplitude close 
to the boundary. 
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Figure 1. Persistent sites (red) in a 1+1-dimensional critical DP process 
(green) . 



Persistence was also studied in various other contexts, including random walks and 
Gaussian processes [9,38-40], fluctuating interfaces and surface growth [41-47], systems 
with shear flow [48-50], models with algebraic long-range interactions [51,52], and 
disordered systems [53-56]. Moreover, various authors have improved and generalized 
the concept of persistence in different ways [57-61,61-68]. Some predictions were 
confirmed experimentally and successfully applied to financial data [64,69-74]. 

Within this broad range of models and applications, the present paper addresses a 
particular issue, namely, the problem of local persistence in the directed percolation 
universality class (DP) [77, 78, 82], the simplest and paradigmatic class of non- 
equilibrium phase transitions from a fluctuating phase into one or several absorbing 
states. In models with absorbing states local persistence may be defined as the 
probability that a given site - monitored shortly after the initial state up to time t 
- remains locally absorbing (inactive). For example, in a ci+l-dimensional DP process 
starting with a fully occupied infinite lattice, Pe{t) would be the density of sites that 
have never been activated up to time t, excluding the initial state at t = 0. A typical 
run of a 1+1-dimensional DP process is shown in Fig. [T], where locally persistent sites 
are represented as vertical red lines. 

In the context of absorbing phase transitions local persistence was studied for the 
first time in Ref. [30], where the fraction of unaffected sites in a DP process was in- 
vestigated numerically. At criticality in 1-1-1 dimensions the persistence probability was 
found to decay algebraically as Pe{t) ~ t~^' with an exponent O^^^ = 1.50(2) which is 
apparently independent of the bulk exponents (3, u^, v^. Later Albano and Munoz [79] 
extended these studies to higher dimensions. Investigating the ZGB model at its DP 
transition point, they reported the estimates Og, = 1.50(1) in one and two, 9i = 1.33(3) 
in three, and Oi = 1.15(5) in four spatial dimensions. 
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In the present work we revisit the problem of local persistence in DP from various 
perspectives. The main results are: 

a) The algebraic decay in Eq. ([T]) is confirmed within a mean field approximation. 
However, in mean field the mechanism leading to a power-law behavior is shown to 
differ from the one in low- dimensional systems. 

b) We provide a more accurate estimate of 6i in 1-1-1 dimension. In higher dimensions, 
however, it is shown that long transients make a reliable numerical estimation 
impossible, questioning previous estimates reported in Ref. [79]. 

c) We discuss the scaling properties of local persistence, suggesting a formulation in 
the continuum limit and investigating the concept of graded persistence. 

d) For the first time local persistence is investigated in seed simulations, where a 
different local persistence exponent is observed. 

e) We analyse spatial (as opposed to temporal) persistence and find a relation between 
the corresponding exponents. 

2. Mean field calculation 

In directed percolation, the mean field approximation is expected to be valid above the 
upper critical dimension dc = 4. However, the results of this section indicate that, as in 
other react ion- diffusion models [5] , the local persistence exponent for dc < d < oo turns 
out to be d-dependent and non-universal. To illustrate this, we shall write down mean 
field approximations for two different lattice models in the DP universality class: the 
contact process [84] and bond directed percolation [85]. 

2.1. Mean field approximation for the contact process 

The contact process is a lattice model with random-sequential dynamics, where particles 
create offspring {A —>■ 2 A) at rate A and self-annihilate {A 0) at rate 1. The mean- 
field evolution equation of this process reads 



where p(t) is the coarse-grained density of active sites. In addition, let ipit) be the density 
of persistent sites that have not been activated up to time t. Since the mean-field limit 
ignores fluctuations, the loss of persistent sites is proportional to the frequency of new 
activations multiplied by the concentration of persistent sites, leading to the additional 
differential equation 



\p={\-l)p{t)-\p\t), 



(2) 



\p{t)^{t) . 



(3) 
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Setting p(0) = ip{0) — 1 the solution at the critical point A = Ac = 1 is given by 

p{t) = m = ^ ^ t-i . (4) 

This means that for the contact process the persistence exponent within mean field 
theory is given by 

er = 1 (5) 

in any space dimension d > dc- It is also instructive to study the mean field solution 
above and below the critical point. Here the solution reads 

In the inactive phase A < 1, the density of active sites decays exponentially while the 
density of persistently inactive sites saturates for large times at a constant: 

p{t) ~ (l-A)exp[-(l-A)t], (7) 
i;{t)^l-X. (8) 

In the active phase A > 1, however, the roles between p and ip are exchanged and one 
obtains 

Pit) - ^ , (9) 

m - ^ exp[-(A - l)t] . (10) 

Since the time scale in the exponential functions is expected to scale as |A — ll"*^!! and 
the density of persistent sites should saturate in the second case at a level proportional 
to |A — ll"^^*^!!, these results are compatible with the well-known mean field exponents 
P^'^ — u^'^ = 1 of directed percolation. 



2.2. Mean field approximation for bond directed percolation 

In order to illustrate in detail that the persistence exponent in finite dimensions above 
dc is generally a d-dependent non-universal exponent, let us now consider the example 
of directed bond percolation. In bond DP on a tilted square lattice each lattice site is 
connected to 2d neighboring sites via directed bonds which are open with probability p 
and closed otherwise. Assuming the states of these neighboring sites to be uncorrelated, 
the probability that the updated site becomes active is 1 — (1 — p)^, where k G {0, . . . 2d} 
is the number of open bonds. Since the probability to find k out of 2d bonds opened 
is given by (^l^)p''{l — p)^'^"^, the average density of active sites after the update, pnew, 
within mean field approximation reads 

^new=f;('^^y(l-pr^-Ml-(l-p)1 

= 1 - (1 - ppf'' = 2dpp - {2d^ - d)p^p^ + 0(p3) . (11) 
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At the mean- field critical point Pc = l/2d this expression expanded to second order 
in the continuum limit reduces to the differential equation p = — 1(1 — ^)p^ with the 
asymptotic solution p{t) ~ — i/2rf)+2 * Gri^^^ this solution cit criticctlity tli6 density of 
persistent sites , to leading order in p, evolves according to the differential equation 

dtm = -lpit)m. (12) 

where the factor | stems from the fact that in bond-DP a given site is updated only 
every second time step. The solution is 



2d-l 

t + 1 



2d 



-2d 



_ 2d 



t^ , (13) 



hence the local persistence exponent is given by 

This result differs from the one for the contact process in Eq. ([5]). It demonstrates that 
above the upper critical dimension, where the mean field approximation is expected to 
give correct answers, 6e is a non-universal mo del- dependent quantity which generally 
depends on the dimensionality. Similar non-universal results were obtained for other 
react ion- diffusion models such as pair annihilation 2A — * 0, see e.g. [5]. 

We finally note that the simple mean field results ([5]), f|T^ rely on the assumption 
that the local density of active sites is spatially homogeneous, i.e., p{x,t) takes the 
same value p(t) both close and far away from persistent sites. In the full model this 
assumption is expected to be valid only in the limit d —>■ oo. In finite dimensions 
above dc, however, the density of active sites next to persistent ones is expected to be 
lower because a persistent site itself can be regarded as a defect of vanishing activity. 
These (anti) correlations may effectively change the density p(t) of active sites near the 
persistent ones, i.e. p{t) = 'jdpit) with 7^ ^ 1 for d ^ 00. Correspondingly, one should 
include this factor in Eqs. ([3]) and ( JT2l) . leading to a further correction to the mean-field 
exponent 9i — > •jd^i- 

The non-universality of 6e for d > dc shows that the local persistence exponent 
differs qualitatively from the other bulk exponents of DP which are constant above 
the upper critical dimension. Apparently the Gaussian fixed point, which governs the 
critical behavior above dc, does not impose constraints on the value of 6£. 



3. Numerical estimation of 6^ 

In Ref. [30] the local persistence exponent of l+l-dimensional DP was determined by 
simulating a l+l-dimensional Domany-Kinzel cellular automaton [82,83], leading to the 
estimate 6e = 1.50(2). This value suggested that Oe is probably unrelated to the other 
bulk exponents of DP, /3, z/y and z/^. Similar estimates were found in other models for 
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Figure 2. High-precision measurement of local persistence in the Id contact 
process. The plot shows P£{t) multiplied by t^^ with 6i = 1.512. The red dotted 
lines indicate the error margin of the estimated exponent ±0.006. The inset 
shows the raw data of the simulation. 




DP [79], suggesting that 9e is universal in 1+1 dimensions. Moreover, it was speculated 
whether 6£ is exactly equal to the rational value 3/2 or not. 

To address this question we performed extensive Monte-Carlo simulations at the 
critical point of the 1+1- dimensional contact process [84], ranging up to t = 10^ time 
steps and averaged over 2.4 x 10^ independent runs. This enormous number of sweeps 
is necessary to obtain reliable averages for the persistence probability which decreases 
down to 10~^. As shown in Fig. [2] the results confirm the anticipated power-law decay 
of Pi{t). Moreover, the local persistence exponent is estimated by 

= 1.512(6) (15) 

so that we can now safely exclude the scenario of a rational value 3/2. A similar 
numerical result was obtained for directed bond percolation with parallel updates. 



model 


bond-DP 


ZGB 


CP 


bond-DP 


ref. 


[30] 


[79] 


this work 


this work 


d = 1 


1.50(2) 


1.50(1) 


1.512(6) 


1.517(10) 


d = 2 




1.50(1) 


> 1.62 


> 1.58 


d = 3 




1.33(3) 


> 1.5 


> 1.4 


d = A 




1.15(5) 




> 1.1 



Table 1. Estimates for in directed percolation. The inequalities for higher 
dimensions rely on the assumption that the curvature of the data shown in 
Fig. [3] does not change sign as t — > oo. 
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The measurement of persistence in higher dimensions is extremely difficuh since 
finite-size effects become increasingly relevant. For example, in three dimensions the 
simulation time of 10^ time steps would require about 10^ lattice sites in order to exclude 
finite-size effects. In addition, a large number of runs is needed to obtain reasonable 
statistical averages. Fig.[3]shows Pe{t) measured in the one-, two- and three-dimensional 
contact process at criticality within a temporal window where finite-size effects can be 
safely excluded. As can be seen, the lines for d = 2,3 bend downward. Therefore, a 
reliable estimation of the critical exponents is impossible; judging from the numerics 
it is not even clear whether Pe{t) exhibits an asymptotic power law at all. However, 
assuming that the asymptotic decay is algebraic and that the curvature does not change 
sign as t ^ oo the simulations allow us to specify lower bounds for 9i, see Table [31 Note 
that the lower bounds for two spatial dimensions are actually larger than the d = 1 
estimate, while they become smaller for = 3, 4, a somehow odd result that may indeed 
indicate abnormally long transient effects. 

To summarize, we obtain a reliable estimate for 6e in 1+1 dimensions while we 
are unable to confirm or disprove a possible asymptotic power-law decay in higher 
dimensions. However, the results rule out a conjectured numerical coincidence of 6e in 
1+1 and 2+1 dimensions (called superuniversality) as previously speculated in Ref. [79]. 

4. Scaling properties 

The observed power-law behavior in 1+1 dimensions suggests to describe the interplay 
of several parameters by phenomenological scaling laws. In this section we discuss such 
scaling laws on different levels of complexity. 
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Figure 4. Local persistence probability in finite-size (left) and off-critical 
simulations (right) in the 1+1-dimensional contact process. The graphs 
demonstrate data collapses according to the scaling form (|16p . while the insets 
show the original data. 



4-1- Off-critical and finite- size scaling 

The simplest extension is to include the lateral system size L and the distance from 
criticality A = A — Ac as additional parameters. In analogy to other quantities exhibiting 
DP scaling one expects an asymptotic scaling law of the form 

Pe{t) ^ r^^F(-^,tA'^ii), (16) 

where F is a scaling function and z = z/||/z/_|_ ~ 1.581 is the dynamical exponent of DP in 
1+1 dimensions. In order to demonstrate the validity of this scaling form numerically, 
we consider two special cases: 

(a) Finite-size scaling: A finite 1+1-dimensional contact process at criticality A = 
reaches the absorbing state within finite time so that there is a finite probability for 
persistent sites to survive forever. Therefore, the persistence probability saturates 
at a constant value. Plotting Pi[t)L^^^ against t/L^ all curves collapse onto a single 
one, as demonstrated in Fig. HI 

(b) Off-critical simulations: For an infinite 1+1-dimensional contact process below 
(above) the percolation threshold, we expect the persistence probability to saturate 
(decay exponentially). According to the scaling form f|T6|) the curves should collapse 
in both cases when Pi{t)/Sr^^^^^ is plotted against tA^^ii. As shown in the right panel 
of Fig. m this is indeed the case. 

The successful data collapses confirm conventional power-law scaling in 1+1 dimensions. 
Whether similar scaling laws apply in higher dimensions is still an open question. 
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Figure 5. Measurement of the persistence probability Pi{t,t()) in a 1+1- 
dimensional critical contact process as a function of two variables t and Iq. 
The left panel shows the original data while the right one demonstrates a data 
collapse according to Eq. (flS]) . 

4-2. Local persistence as a quantity depending on two time parameters 

In models with absorbing states local persistence is usually introduced as the probability 
that a given site is not activated until time t, excluding the initial configuration at 
t = where all sites are active. This definition implicitly suggests that the persistence 
probability Pe{t) depends on only a single time parameter t. However, with only one 
parameter t it is impossible to define a meaningful continuum limit of the persistence 
probability. For example, a fully occupied initial state in the lattice model translates 
into an infinite 5-peak-like density of active sites at t = in the continuum description 
so that the probability for a given site to remain inactive in the interval < r < t would 
be zero. Therefore, the only way to introduce persistence consistently is to define Pe as 
the probability that a given site remains inactive in some interval between time to > 
and time t. This means that local persistence should be defined as a quantity Pi(t,to) 
that depends on two time parameters to and t which reduces to the previous definition 
by setting to := 1- 

The left panel of Fig. [5] shows the temporal decay of the persistence probability 
Pe{t, to) in a critical 1+1-dimensional contact process as a function of t for various values 
of to. As can be seen, by increasing to the initial transient becomes more pronounced. 
Moreover, the curves are shifted to the right. However, the asymptotic slope, i.e., the 
persistence exponent 9e, is not affected. 

In order to describe this behavior by appropriate scaling laws, the straight-forward 
method would be to include both parameters to and t in the scaling form and to require 
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invariance under scale transformations of the form 

A^AA, p^A^p, L^A-'^^L (17) 

to ^ A-'^ii to , A-'^ii t, ^ A"P,, 

where A > is a scale factor, A = A — Ac is the distance from criticality, p denotes 
the density of active sites, L is the system size, and k is an exponent to be determined 
below. At criticality in an infinite system, invariance under this set of transformations 
would lead to the scaling form P^(t, to) = t^'^/^'H P(t/to) with an unknown exponent k. 
This exponent can be determined by keeping the difference At = t — to fixed and taking t 
to infinity. In this limit the critical system approaches the absorbing state and therefore 
the probability Pi(t,to) to find no activity between to and to + At tends to 1. This 
implies that /t = 0, i.e., the local persistence probability is a quantity that does not 
carry an intrinsic scaling dimension. Therefore, we arrive at the scaling form 

P,(t,to) = P(t/to), (18) 
where P is a scaling function with the asymptotic behavior 

F{z) ~ z-^' for z > 1, (19) 
F{z) ^1 for 2 ^ 1. 

It is interesting to note that autocorrelation and autoresponse functions in non- 
stationary critical systems behave similarly. These two-time quantities were studied 
in recent years in the context of aging [86] and exhibit similar scaling laws, of course 
with different exponents and scaling functions. 

In order to verify the scaling form ( ITSl) and its asymptotic behavior (|T9l) we 
simulated a 1+1-dimensional contact process and measured the probability that a given 
site is not activated in the time interval between to and t. The left graph in Fig. [5] shows 
curves for fixed to as a function of t. The graph on the right side demonstrates the data 
collapse according to the scaling form (fTSjl . plotting Pi{t, to) versus the quotient t/to — 1. 
Although the collapse seems to be convincing, it is not really perfect. The reason will 
become clear in the following subsection. 



4-3. Graded persistence 

In 1996 Ben-Naim et al. [88] introduced a parametrized generalization of local 
persistence, namely, the probability Pn(t) that a given site changes its state exactly 
n times until time t. Clearly, for n = one retrieves the usual definition of local 
persistence. Studying the voter model the authors showed that this quantity obeys the 
scaling form 

p„w^ijt(^)- 

Later, Majumdar and Bray [59] considered a similar generalization called persistence 
with partial survival P{p,t), where upon a spin fiip a site is decleared as non-persistent 
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with probability p. This definition reduces to the original one in the limit p ^ 1 and is 
related to the previous generalization by P{p,t) = Yl'^=oP"'^n{t)- 

In addition, Majumdar and Bray argued that there are in principle two different 
scenarios: If the underlying process is 'smooth', meaning that the average number of zero 
crossings (spin fiips) per site (n) is finite, one has P{p, t) ~ t~^^ with a generalized local 
persistence exponent 6p that varies continuously with p. This situation is observed e.g. 
in the diffusion equation starting with a random intial field. On the other hand, if the 
underlying process is 'non-smooth', meaning that (n) is infinite, there is no continuously 
varying persistence exponent, instead the scaling behavior discussed by Ben-Naim et al. 
is recovered. This happens e.g. in the Glauber-Ising model. In the case of directed 
percolation the average number of spin fiips per site is directly linked to the integrated 
particle density {n) ~ dt p(t), hence {n) is infinite below the upper critical dimension 
so that the scenario of a 'non-smooth' process should apply. 

It is interesting to note that for DP such a generalized persistence probability can 
also be motivated by field-theoretical arguments as follows. In a continuum formulation, 
the discrete variables of the lattice Si{t) are replaced by a coarse-grained density p(r,t) 
of active sites. As suggested by Cardy [5], a suitable definition of the local persistence 
probability in the continuum limit is given by 

P,{t,to,L,A,fi) = ^exp(^-i ^ dt'p{f,t')^^ . (21) 

In this expression, which is also needed to describe epidemic process with 
immunization [89-94] , an additional parameter p has to be introduced in order to make 
the argument of the exponential function dimensionless. This parameter plays a similar 
role as n in Eq. fl20|) and may be interpreted as a typical threshold for the integrated 
activity below which the site at position r is 'declared' as being persistent. 

Invariance of the argument of the exponential function under scale transforma- 
tion (fT71) requires that the parameter p scales as 

p ^ A^-^Wp (22) 

so that Pi is a generalized homogeneous function of the form 

Pe{t, to, L, A, p) = PiiA-^H, A-'^Ho, A'^^^L, AA, A^^^H/i) . (23) 

Therefore, a scaling form for the local persistence probability, similar to the one in 
Eq. fl2U]) . is given by 

Pi{t,to,L,A,p) = f(-, ^, -V, ^7^^) , (24) 

where a = (3/v\\. For L oo, A ^ 0, and p ^ this scaling form reduces to 
Eq. ffTSjl . which would predict a perfect data collapse in Fig. [51 However, in most 
numerical simulations, p is effectively given by a small but finite number of order 1. 
Thus, dropping the last argument in Eq. will lead to small errors. This is probably 
the reason why the data collapse shown in Fig. \5\is quite satisfactory but not perfect. 




Figure 6. Graded persistence. The figure shows measurements of the 
probability (t) that a given site is activated less than m times at time t. 
The inset shows the original data while the main graph demonstrates the 
corresponding data collapse. 



In order to verify this scaling law quantitatively, we consider a graded persistence 
probability 

P/"^)(t,to) = (x(m-^..(t'))) , (25) 

\ t'=to I 

where Siit) = 0, 1 is the state of site i and 

^«^{; 

is a step function. Obviously P/"^\t, to) is the probability that a given site i is activated 
less than m times during the time interval from to until t. For = it is related to 
the generalization by Ben-Naim et al. by P/"^\t, 0) = Yl^=o Pn{t)- The reason why we 
define graded persistence as the probability for having less than m activations is that 
the parameter m G N plays the same role as the constant /i in Eq. (pT]) . 

To verify the postulated scaling behavior in Eq. fl24l) we measured Pi^'^\t,to) for 
to = 1 in a critical system for various values of m. As demonstrated in Fig. El by plotting 
Pi^"^\t) versus t/m}/^^~^'^ all curves collapse onto a single one. This data collapse 
confirms a posteriori the field-theoretic definition of local persistence in Eq. (I2TI) . 
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Figure 7. Rapidity reversal symmetry: Reversing the time arrow the local 
persistence probability Pe{t) in bond-DP is exactly equal to the return 
probability of a cluster generated by a permanently active seed (see text). 

5. Mapping local persistence to the problem of a permanently active site 

The field-theoretic action describing DP transitions is known to be invariant under the 
so-called rapidity reversal symmetry [87], implying that the density of active sites p{t) 
in systems with homogeneous initial conditions and the survival probability P{t) in seed 
simulations exhibit the same type of asymptotic power-law scaling. In bond directed 
percolation, this symmetry even holds exactly on the microscopic level since for a given 
realization of open and closed bonds the probability to find a directed path from site A 
to site B forward in time and the corresponding probability to find a directed path from 
B to A backward in time are identical. In the following let us consider the implication 
of the rapidity reversal symmetry to the problem of local persistence in DP. 

5.1. Relation to a DP process with a permanently active source 

As illustrated in left panel of Fig. U\ a site is locally persistent (vertical red line) if it 
is not reached by any directed path from the initial state (upper horizontal green line) 
until time t. This means that the cluster generated by the initial state (green pixels) 
and the set of all directed paths leading to the persistent site (red pixels) are disjoint. 

The right panel of Fig. [7] shows exactly the same realization of open and closed 
bonds turned upside-down which can be interpreted as a bond-DP process running 
backward in time. As can be seen, the red cluster may now be interpreted as being 
generated by a permanently active site in the middle. Obviously, the local persistence 
probability is exactly equal to the probability R{t) that this cluster does not reach any 
site at the final time (horizontal blue line). In other words, Pe{t) coincides with the 
return probability R{t) of a DP cluster generated by a permanently active site in an 
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otherwise empty system to its initial condition (i.e. an empty lattice but the single 
permanently active site) after t time steps. 

In other realizations of DP, e.g. in the contact process, the time reversal symmetry 
does not hold exactly on the microscopic level, instead it emerges as an asymptotic 
symmetry on large scales, just in the same sense as a random walk on a square 
lattice becomes rotationally invariant on large scales. Therefore, apart from a non- 
universal proportionality constant, the relation between local persistence and the return 
probability of clusters generated by a permanently active seed is expected to be valid 
in all models belonging to the DP universality class [30]. 

5.2. One-sided persistence in one spatial dimensions 

In one spatial dimension, the permanently active seed splits the lattice in two disjoint 
uncorrelated domains, so that the return probability to the initial condition in the 
entire space is equal to the square of the one-sided return probability when only a single 
semispace is considered, i.e. Rh{t) ~ t"^*^/^. In principle, this would allow for a better 
precision in numerical simulations, since a slower decay would lead to a better statistics 
for large times. A measurement of Rh{t) at the critical point of bond directed percolation 
is shown in Fig. [HK, confirming the halved value of the persistence exponent. In higher 
spatial dimensions, the permanently active seed no longer splits the space into disjoint 
domains, and so one has to deal with the return probability to the initial condition in 
the full space. 

5.3. Return vs. first-return probability 

In order to verify the identity between local persistence probability and the probability 
to return to the initial condition of clusters generated by a permanently active seed 
in higher dimensions, we chose to measure the first-Tetuin probability distribution Qi 
which is defined as follows: Consider an active seed located at the origin of a lattice in 
d spatial dimensions. Let Qi(t) be the probability for the cluster generated from the 
active seed located to return to its initial condition at time t for the first time, with 
Qi{t) = for t <= 0. Then it follows that the probability Q2(t) of the second return to 
the initial state, the last one at time t, is given by the discrete convolution 

oo 

Q2(t) = J]QiHgi(t -n) = (Qi * Q,){t) (27) 

n=l 

since successive returns to the initial state are uncorrelated. This can be extended 
recursively to the probability of k returns (with k > 2), i.e. Qk(t) = {Qk-i * Qi)if)- 
Suppose that Qi{t) ~ t~'^ with 7 > 1. By the convolution theorem for the Laplace 
transform one has the asymptotic Qk{t) ~ It is now easy to verify that the 

unrestricted return probability R{t) can be expressed as a sum of the k first-return 
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Figure 8. (a) One-sided return probability Rhit) of clusters generated by a 
permanently active seed for critical bond directed percolation confined in the 
semiplane of a (l-l-l)-dimensional tilted square lattice (see text). The dashed 
red line marks the power low decay t"^*?/^ with 6g = 1.517. Data has been 
averaged over 4 • 10^ independent realizations, (b) First-return probability 
distribution Qi{t) of clusters generated by a permanently active seed for critical 
bond directed percolation on a (2-|-l)-dimensional body-centered cubic lattice. 
The red dashed line indicates the lower bound for 6^ in two spatial dimensions 
reported in Table [3l Data has been averaged over 2.5 • 10^ independent 
realizations. 



probability distributions, 

t 

R{t) = Y,Qk{t) (28) 

k=l 

which scales asymptotically as t~'^. Therefore, to leading order, from the rapidity 
reversal argument one simply expects 

l = 9e. (29) 

In Fig. [H)d the decay for the first-return probability Qi(t) at the critical point of two 
spatial dimension bond directed percolation is compared with the lower bound reported 
in Table [31 



6. Local persistence in seed simulations 

Up to now we considered only a fully occupied lattice po = 1 as initial configuration. 
More generally one could start with a finite density < po < 1 of active sites. In [30] 
it is shown that in 1-|-1- dimensional DP, after an initial transient, the local persistence 
probability follows the power law P^it) ~ t^^* with an exponent 6^ which is independent 



Local Persistence in the Directed Percolation Universality Class 



17 




Figure 9. Local persistence in seed simulations. Left: Lattice configuration 
in (1+1) dimensions. Each lattice site has two nearest neighbors, which could 
be activated in the next time step, following the red path. A first-return to 
the origin is only possible at even times, starting with one single seed at the 
origin at time t = (green). Middle: Snapshot of the evolution of a surviving 
percolation cluster in one spatial dimension. Only the right half of the cluster 
is active. Right: Lattice configuration in (2+1) dimensions. Again, possible 
paths (red) and the seed at the origin (green) are shown. 

of Pq. In the limit po —>■ the transient regime takes an infinite amount of time, meaning 
that asymptotic power law is no longer accessible. 

Another well-known simulation technique starts with a single 'seed' at the origin of 
the lattice, which generates a certain cluster of active sites (see Fig. IH]). This simulation 
technique was first introduced in 1979 by Grassberger and de la Torre [87]. Many of these 
clusters will die out soon, but some of them will survive for long time. The advantage of 
this approach is that the number of particles at each time step is finite. Thus, instead 
of storing the whole lattice configuration one only has to keep track of the coordinates 
of the active sites in a dynamically generated list. Apart from the range of integer 
numbers, this approach imposes no boundary effects, thus one effectively simulates an 
infinite system. 

It is interesting to investigate how the local persistence probability Pl'""^{t) behaves 
in simulations starting with a single seed. In this case local persistence may be defined 
as the probability that the origin, from where the cluster is generated, is not activated 
until time t. However, since a finite fraction of clusters dies before they have returned 
to the origin, this quantity, when averaged over all runs, would tend to a constant for 
t — s> oo. Therefore, local persistence should be defined as an average over the surviving 
clusters. The question is, whether 

(i) P['""^{t) follows a power-law behavior as Pi{t) does in a system with homogeneous 
initial conditions, and 

(ii) if so, whether and how the corresponding exponent 6/"'"^ is related to 6e. 

Equivalently one can measure the probability distributioiiji Df'"^{t) for the first return 
I In [79], which deals with homogeneous initial conditions, the corresponding quantity Di{t) is called 
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Figure 10. Local persistence in seed simulations. Left: Time distribution of the 
frequency N{t) of the first-return in 1 (green), 2 (red), 3 (blue) and 4 (yellow) spatial 
dimensions. For better visualization neighboring data points are joined. Right: One- 
dimensional data Df'"^{t) divided by the fitted power law Dl'''°'"^{t) cx t^^' plotted 
against time t (green). The red and the blue line indicate the error margins for the 
estimate of 9g'"""^. 



Note that Df^^lt) differs from the first return probabihty distribution Qi{t) considered in 
the Sect. [5] in so far as the latter refers to the return to the original initial condition when 
the seed is kept permanently active in time, while the former samples first-reactivation 
times of a normal seed. 

Assuming a power-law behavior, we expect a different value for the local persistence 
exponent ^^"'""^for the following reason: The initial condition of a single seed can be 
interpreted as a homogeneous initial state in the limit po 0, with an important 
difference: We only investigate the first-return probability of a single seed at the origin, 
which is a non-homogeneous problem. In a sense, we exclude all lattice sites except the 
origin from the calculation of the first-return probability. Consequently, Df^^lt) < Di{t) 
(for all t), which immediately leads to Oi""""^ > 6e. In addition the clusters have a finite 
lifetime so that Df'^lt) is again expected to decline faster than Di{t). Therefore it is 
natural to expect Oi""""^ to be greater than 6e. Moreover, persistence in seed simulations 
is not defined in the limit of infinitely many dimensions since then the generated cluster 
never returns back to the origin. 

In order to verify this conjecture we simulated the bond-DP process on a simple 
cubic lattice in c? = 1 to 4 spatial dimensions. In this model, each lattice site has 2d 
nearest neighbors. We start with a single seed at the origin. Notice that because of 
the lattice structure, a first-return to the origin can only occur at even time steps (see 
Fig. [9]). Many of the evolving percolation clusters die out before a first-return to the 
origin occurs (approximately 49% in one, 56% in two, 74% in three and 83% in four 
spatial dimensions). We simulated about 10^° independent runs in dimensions two to 




(30) 



the persistence probability distribution, which gives the probability for a persistent site to become active 
in the time interval between t and t + dt. 
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four and, in order to increase the precision of er^"'''', about 10^^ independent runs m 
the one-dimensional case. 

The results are displayed in Fig. [TOl In the one- dimensional case the evaluation 
of the data yields a clean power-law for the first-return-probability distribution 
with exponent 6^^'""^'^'^ + 1 = 3.014±0.012. Although the integer value 3 lies only slightly 
outside of the error margins and we therefore cannot exclude this value for sure, we see 
no reason for Q^'^'^'^'^^^^ to be integer. The higher- dimensional cases also appear to follow 
a power-law behavior. However, a more detailed analysis by plotting the data divided 
by the fitted power-law shows systematic deviations. If we anticipate a power-law, limits 
for the corresponding exponents are given by 0/"'"^'^° + 1 > 2.5, + 1 > 2.2 and 

^^.eed,4D ^ ^ ^ 2. The questlou whether Or"" is related to Oe is still open. 

7. Temporal vs. spatial persistence 

So far we defined local persistence as a time-dependent probability Pe(t) that a given 
site remains inactive from time to = 1 until time t. For homogeneous initial conditions 
Pi{t) decays as t~^% where 9^ is the persistence exponents. Since space and time play a 
similar (although not a symmetric) role in DP, the question arises whether it is possible 
to consider a spatial persistence probability P(r) with analogous properties. Obviously, 
such a spatial persistence probability should be defined as the probability that a line of 
sites extending from the origin to site r next to an active seed is found to be inactive. 
The problem of spatial persistence was first introduced and analyzed by Majumdar and 
Bray in the context of fluctuating interfaces [95]. 

As shown in Fig. [11] spatial persistence in DP is most easily introduced in 1+1 
dimensions. Part (a) of the figure shows the one-sided version of ordinary temporal 
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persistence: Starting with an initial state where all sites to the right of the origin are 
occupied, the one-sided persistence probability Ph{t) is the probability that no path 
from the horizontal line of initially occupied sites reaches the vertical line at the origin 
extending up to time t. This probability is known to decay as Phit) ~ where 
is the local persistence exponent. 

Part (b) shows the same arrangement rotated by 90°: Here we ask for the probability 
P/j(r) that there is no path from a vertical line at the origin — oo < t < to the finite 
horizontal line extending from the origin to site r. More specifically, we keep the site at 
the origin active and let the system evolve. Once it reaches a stationary state we ask 
for the probability that the sites 1,2, ... ,r are inactive. Obviously, this quantity is a 
spatial analogue of temporal persistence. 

Fig. [12] shows a numerical measurement of i\(r) in the stationary state of a critical 
contact process with an active seed at the origin. We observe that the one-sided spatial 
persistence probability decays algebraically as 

Phir) ~ r-""/^ , k/2 = 1.205(8) (31) 

This estimate is consistent with the scaling relation 

K = eez, (32) 

where z = i^w/vi. is the dynamical exponent of DP. Clearly, the two-sided spatial 
persistence probability that r sites on both sides of the active seed are inactive will 
scale as P(r) ~ r~'^ . 

In dimensions d > 1 spatial persistence may be defined as the probability in a 
stationary DP process with an active seed at the origin that all sites in a c?-dimensional 
sphere of radius r surrounding the origin are inactive. 

8. Concluding remarks 

In this work we have analyzed in detail different aspects of local persistence in DP 
processes. Concerning persistence in d > = 4 spatial dimensions, a simple mean field 
approach shows that the local persistence exponent 9i is a non-universal quantity related 
to the amplitude of the density field, which generally depends on both dimensionality 
and model details. A similar non-universal behavior above the upper critical dimension 
was previously observed in other reaction diffusion models as well [5]. 

Below the upper critical dimension, we have performed high precision numerical 
simulations in order to improve current estimates of 9i in one spatial dimension. On this 
basis we can safely exclude the conjecture 6^ = 3/2, i.e., the local persistence exponent 
- like other DP critical exponents - is seemingly not given by a simple rational value. 

The simulations in c? = 2,3,4 spatial dimensions are less conclusive, apparently 
due to huge transient effects, which make it impossible to estimate the value of 6i 
by numerical simulations. However, specifying lower bounds we can at least rule out a 




Figure 12. Monte Carlo simulation of spatial persistence in a 1+1-dimensional contact 
process at criticality. 

possible superuniversality between d = 1 and d = 2 mentioned in [79]. Although we have 
not been able to determine Of in d < 4, or even to prove the existence of an asymptotic 
power law behavior analytically, we have discussed various concepts - scaling relations, 
graded persistence, mapping to active seed problems and spatial persistence - which we 
hope could prove useful for further theoretical and numerical advances. 

Finally, let us add a comment on our results concerning local persistence in 
seed simulations. Numerical results exhibit a clean power-law in the one-dimensional 
case, whereas in higher dimensions we observe systematic deviations from a power-law 
behavior within the numerically accessible range of time. It is not clear whether these 
deviations just resemble a longer-lasting transient period. They could also be related to a 
qualitative difference between the one-dimensional and the higher-dimensional cases, or 
more precisely to a kind of special screening effect in one dimension: In one dimensional 
seed simulations the persistence of the lattice site at the origin is not influenced by 
remote branches of the emerging cluster, because these are shielded by the branches next 
to the origin (left and right), as long as those have not died out. If a remote branch 
approaches the origin, it will merge with the inner branches without changing their 
dynamics (this is particularly easy to see in l-|-l-dimensional directed bond percolation). 
In other words, a return of the seed to the origin can only be accomplished from one 
of the two branches next to the origin. In higher dimensions, there is no such shielding 
effect. Take for example the (2+l)-dimensional case: The initial seed develops into a 
cone-like cluster, and the origin could be reached from any direction in the plane of the 
lattice, not only from the nearest branches. So, with the breakdown of the shielding 
effect, a return to the origin is more likely, which leads to a stronger, perhaps even 
non-algebraic dechne of local persistence. 
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